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Abstract 

Effective field theory predicts that the leading nonuniversal effects in the 
homogeneous Bose gas arise from the effective range for S-wave scattering 
and from an effective three-body contact interaction. We calculate the lead- 
ing nonuniversal contributions to the energy density and condensate fraction 
and compare the predictions with results from diffusion Monte Carlo calcula- 
tions by Giorgini, Boronat, and Casulleras. We give a crude determination of 
the strength of the three-body contact interaction for various model poten- 
tials. Accurate determinations could be obtained from diffusion Monte Carlo 
calculations of the energy density with higher statistics. 
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I. INTRODUCTION 



The Bose-Einstein condensation of trapped atoms allows the experimental study of cold 
Bose gases with high precision. It is well-known that the dominant effect of the interactions 
between the atoms can be captured by a single constant a called the S-wave scattering length. 
This property is often called universality. However, sufficiently accurate measurements will 
reveal the sensitivity to aspects of the interatomic interactions other than the scattering 
length. These are called nonuniversal effects. 

Theoretical investigations of the homogeneous Bose gas in the 1950's showed that its 
properties could be calculated using a low-density expansion in powers of V na 3 , where n is 
the number density. For example, the energy per particle has the expansion 
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The leading term is the universal mean- field contribution first determined by Bogoliubov [|T| . 
The V na 3 correction was first calculated by Lee, Huang, and Yang for bosons interacting 
through a hard-sphere potential 0. In the na 3 correction, the coefficient of log(na 3 ) was first 
calculated by Wu for a hard-sphere potential pi. These leading terms in the expansion were 
shown to be universal applying to bosons interacting through any short-range potential 
with scattering length a. On the other hand, Hugenholtz and Pines showed that the constant 
d in (TJ) as well as higher-order corrections are not universal [[| . They depend on properties 
of the interactions between the bosons other than the S-wave scattering length, although 
the specific properties were not identified. 

Giorgini, Boronat, and Casulleras have studied the properties of the ground state of 
homogeneous Bose gases numerically using a diffusion Monte Carlo method |J. They con- 
sidered 4 different two-body potentials as models for the interaction between the bosons, and 
calculated the energy per particle and condensate fraction for the homogeneous gas. In the 
case of the energy per particle, the universal Vna 3 correction in (jl|) accounts for most of the 
deviations from the universal mean-field prediction at the densities studied, which ranged 
from na 3 = to na 3 = 0.244. However, they also observed small differences between the 
model potentials, i.e. nonuniversal effects. 

Effective field theory provides a powerful method for analyzing the nonuniversal effects in 
the homogeneous Bose gas. It identifies precisely which aspects of the interactions between 
low-energy bosons enter at each order in the expansion in Vna 3 . For example, the coefficient 
d in the na 3 correction in (JJ) is determined by an effective three-body contact interaction 
between low-energy bosons [[!]]. The (na 3 ) 3 / 2 correction is completely determined by d and 
by the effective range for S-wave scattering. 

The nonuniversal corrections in the low-density expansion of the Fermi gas at zero tem- 
perature have been studied previously |8HiT|. The low-density expansion can be expressed 
as an expansion in kp, where kp is the Fermi wavenumber. For fermions with a single spin 
state, the leading correction to the energy of an ideal Fermi gas is proportional to kpa 3 , 
where a p is the P-wave scattering length || . For fermions with g > 1 degenerate spin states, 
the universal corrections to the energy can be expanded in powers of kpa, where a is the 
S-wave scattering length ||. The leading nonuniversal corrections are the kpa 3 term and a 



2 



k F a 2 r s term, where r s is the effective range for S-wave scattering. At order kpa 4 , there is 
also a nonuniversal correction from the effective three-body contact interaction ||. These 
corrections were recently discussed within the effective field theory framework by Hammer 
and Furnstahl |TT] . 

In this paper, we use effective field theory to calculate the leading and next-to-leading 
nonuniversal corrections to the energy density and condensate fraction for a homogeneous 
Bose gas. We compare the calculations with the diffusion Monte Carlo data from Ref. |6| 
and attempt to determine the coefficient of the three-body contact interaction for each of 
the interaction potentials that were considered. 

In Section [H], we describe the general philosophy of effective theories for dealing with 
the low-energy behavior of physical systems. We use an effective quantum field theory to 
identify those parameters that give the leading nonuniversal contributions to the properties 
of a homogeneous Bose gas. The calculations of the leading and next-to-leading nonuniversal 
contributions to the energy density and the condensate fraction are presented in Section |TTT]. 
In Section [TV], we describe the model interaction potentials used by Giorgini, Boronat, and 
Casulleras and review their diffusion Monte Carlo results. In Section [V], we analyze the 
diffusion Monte Carlo results of Ref. and give a crude determination of the three-body 
contact parameter for each of the potentials that was considered. A summary and our 
conclusions are presented in Section [VI[ 



II. EFFECTIVE FIELD THEORY 

Effective theory is a general approach to understanding the low energy behavior of a 
physical system that has deep roots in several branches of physics. Some of these roots are 
described in Ken Wilson's Nobel lecture on the Renormalization Group fl2 |. In elementary 



particle physics, the roots have two main branches. One branch evolved from efforts to 
make intuitive sense of renormalization in quantum electrodynamics (QED) fl3| . In per- 
turbative calculations in QED, intermediate steps are plagued by "ultraviolet divergences" 
that indicate strong sensitivity to physics at extremely short distances. The renormalization 
procedure eliminates the ultraviolet divergences and allows extremely accurate predictions 
of the properties of electrons and their bound states in terms of two fundamental param- 
eters: the electron mass and the fine structure constant. The other main branch evolved 
from efforts to understand the low-energy behavior of pions and nucleons. The fundamental 
description of these strongly interacting particles is provided by quantum chromodynamics 
(QCD). However, an effective theory can be used to describe the low-energy behavior of 
these particles accurately without using any information about QCD other than its global 
symmetries and measurements of a few low-energy observables Jl4| . Starting from these two 



main roots, effective field theory has developed into a universal language for modern elemen- 
tary particle physics |TIJ. It is used to develop low-energy approximations to the Standard 



Model of elementary particle physics and also to quantify corrections to the Standard Model 
arising from a more fundamental theory, such as a unified field theory or string theory. 

Most of the applications of effective theories to date have been carried out within the 
context of quantum field theory. However, the principles of effective theory apply equally 
well to the simple quantum mechanics problem of two particles interacting through a central 
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potential U(r) [[Uj. Suppose we are interested only in the low-energy observables of the sys- 
tem, where "low-energy" refers to energy E close to the threshold E = 0. These observables 
include bound-state energy levels close to threshold and the low-energy scattering ampli- 
tudes. Suppose also that we know the potential U(r) accurately at long distances r > R, 
but that its behavior at short distances r < R is not known accurately (If it is a short-range 
potential with range smaller than R, then U(r) = for r > R. If the particles are real atoms 
interacting at long distances through a van der Waals potential, then U(r) = —Cq/t 6 for 
r > R.) Given sufficiently precise information about some low energy observables, effective 
theories allow other low-energy observables to be calculated with arbitrarily high accuracy 
without having any information about the short-distance potential. 

The method is very simple. Simply replace U(r) by an effective potential U e s{r) that 
is identical for r > R and whose form for r < R involves an adjustable parameter c\. 
For r < R, the effective potential need not bear any resemblance to the original potential 
U(r) as long as it has an adjustable parameter. A sufficient condition on the short- distance 
potential is that by varying c\, one should be able to change the number of bound states 
by ±1. This guarantees that the scattering length a can be tuned to any value from — oo 
to +oo. A simple example is a constant potential with adjustable depth V\. Tune the value 
of this parameter c\ so that the scattering amplitude at threshold is reproduced exactly. 
Then the Schrodinger equation with U e s(r) will reproduce all other low-energy observables 
with errors that scale linearly with E. To achieve higher accuracy, use an effective potential 
U e s(r) with two independent adjustable parameters c\ and c-i- A simple example is a 2-step 
potential with adjustable depths V\ and Vi- The equality of U e s(r) and U(r) for r > R 
guarantees that the difference between their scattering amplitudes is an analytic function of 
the wavevector k. If we tune c\ and c 2 so that the S-wave scattering amplitudes agree to 
second order in k = |k|, then U e g(r) will reproduce all low-energy S-wave observables with 
errors that scale like E 2 . By tuning a third parameter C3, we can guarantee that P-wave 
observables also have errors that scale like E 2 . If we add more adjustable parameters and 
tune the scattering amplitude to higher order in E, the errors will scale as higher powers 
of E. Low-energy observables at energies within the appropriate radius of convergence can 
therefore be calculated to any specified accuracy by adjusting a finite number of parameters 
in the effective potential. 

The problem of a Bose gas at zero temperature is an ideal application of effective field 
theory. Suppose we have a gas with number density n consisting of identical bosons of mass 
m interacting through a short-range two-body potential U(r). (For convenience, we will 
refer to these particles as "atoms.") Using the method of second quantization, this system 
can be described by a quantum field theory. The quantum field operator ip(r,t) annihilates 
an atom at the point r, and i/j'i/j is the number density operator. The dynamics of this 
quantum field theory can be summarized by a nonlocal Lagrangian density: 
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- - —Vip* ■ Vip + fM^i/j - - I d 3 r'i/j*ijj(r)U(r - r')f^(r') . (2) 
2m 2 



The parameter [i is the chemical potential. If we wish to describe a system consisting of a 
small number of atoms, such as a two-body or three-body system, then the ^ip*ip term can 
be ignored. To describe a homogeneous gas of number density n, \i must be tuned so that 
the expectation value of iftty in the ground state of the quantum field theory is equal to n. 
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Now suppose we are interested only in low-energy observables, such as the properties of 
the ground state. The principles of effective theory tell us that the short-range potential U(r) 
can be replaced by any other set of effective short-range interactions, so long as they have 
adjustable parameters that can be tuned to reproduce a set of low-energy measurements. 
An example of such a set of low-energy measurements is the coefficients of the expansions 
in y/na 3 for the energy density and other properties of the many-body system. A simpler 
example is the coefficients in the low-momentum expansions for the scattering amplitudes 
of atoms in the two-body sector, three-body sector, etc.. 

If we wish to carry out analytic calculations, one particularly convenient choice for the 
effective short-range interactions is to take them all to be purely local. The dynamics of the 
resulting quantum field theory is then described by a local Lagrangian density: 

1 r ■ -i h 2 

1 i/j*tp Vtp* ■ Vip + u^ip (3) 

J 2m 



-^(V^) 2 - ^v(^) • v (v>v) - §(^) 3 + • • • 

The terms in the Lagrangian are constrained by the symmetries of the original Lagrangian: 
Galilean invariance, parity, and time reversal. We have included explicitly those terms that 
contribute the leading nonuniversal corrections to the properties of the homogeneous gas. 
All other terms are either higher order in ip or in V. The [V (?/>*■?/>)] 2 term is the only possible 
term that is fourth order in ip and second order in V. Up to total divergences, the only 
other independent term is (ip*Vip — ipVip*) 2 , but it is forbidden by Galilean invariance. 
Note that while the original Lagrangian (fj) included only terms of 2nd and 4th order in if), 
the effective Lagrangian (fj) includes terms of 6th order and higher. The 6th order terms 
are necessary, because a local Lagrangian without such terms could not reproduce the low- 
energy behavior of scattering amplitudes in the three-body sector. It would also be unable 
to reproduce the na 3 correction term in the low-density expansion ([!]) for the energy density 
of the homogeneous gas. 

The coefficients Aq and A 2 in @ can be determined by demanding that the effective 
theory reproduces the amplitude for low-energy atom-atom scattering up to errors that 
scale like the square of the energy. The T-matrix element at tree-level for two atoms with 
wavenumbers k 1; k 2 to scatter into atoms with wavenumbers k'^k^ can be deduced from 
the interaction terms in (|3]): 



(T 2 ) tree = -A + -A 2 [(kx - k£) ■ (k 2 - k' 2 ) + (kx - k 2 ) ■ (k 2 - ki)] • (4) 



Going to the center of momentum frame with k x = — k 2 = k and k' x = — k' 2 = k' and using 
the energy conservation condition k 2 = k' 2 , this reduces to 

(7" 2 ) trcc = -A - 2A 2 k 2 . (5) 

The corrections from loop diagrams depend on the regularization and renormalization pre- 
scriptions, but they form a geometric series that can be summed analytically. If we impose 
an ultraviolet cutoff A on the wavenumber, the exact T-matrix element can be determined 
analytically by solving a simple integral equation and has the form 
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where Ao and A2 are complicated functions of A , A 2 , and A. The absence of any angular 
dependence indicates that there is S-wave scattering only. We get a much simpler expression 
if we use dimensional regularization, which automatically subtracts ultraviolet divergences 
that grow like a power of A. In this case the exact T-matrix element is 
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The expressions @ and (|7|) for T 2 do not have precisely the same functional depen- 
dence on k. However, if we make the identifications 1/Aq = 1/Aq + mA / (Airh 2 ) and 
A 2 dT ^ / (A^) 2 = -A 2 /(Ao) 2 , the expansions of © and (0) for small k are identical through 
order k 3 . The difference scales like E 2 , where E = h 2 k 2 /m is the total energy. From the 
effective theory point of view, including the A 2 term in the effective Lagrangian (||) allows 
low-energy observables in the two-body sector to be reproduced up to errors that scale like 
E 2 . Different regularizations require different values of Aq and A 2 , but they give the same 
low-energy observables up to errors that scale like E 2 . 

The amplitude for atom-atom scattering in the original theory described by the La- 
grangian (Q) can be determined by solving the Schrodinger equation with the appropriate 
interatomic potential U(r — r'). The contribution from S-wave scattering can be expressed 
in a manifestly unitary form in terms of the S-wave phase shift 5o(k): 
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The scattering length a and the effective range r s for S-wave scattering are defined by the 
low-momentum expansion of the phase shift. Because of the identity e lS sin 5=1/ (cot S — i), 
the low-momentum expansion is conveniently written in the form 



k cot [d>o(fc)] 
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which defines the effective range r s . After inserting this into @, the T-matrix element can 
be expanded in powers of k. The expansions of (0) and @ match through third order in k 
if the coefficients A^ and A^ have the values 
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These are the coefficients if we use dimensional regularization to remove ultraviolet diver- 
gences. If we used a wavenumber cutoff as in (|), the coefficients A and A 2 in the La- 
grangian would depend on A in a very complicated way. However, their values are uniquely 
determined by the dimensional regularization coefficients ( |T0"D and ([TT|). 
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The (ip*ijj) 3 term in the Lagrangian (|3]) is necessary in order for the effective theory to 
reproduce the amplitude for the low-energy scattering of three atoms with errors that scale 
linearly in the energy E. At tree level, this term gives a momentum independent contribution 
— B Q to the T-matrix element T 3 for 3^3 scattering. One can see the necessity for such 
a term by considering the perturbative expansion for 7-j in powers of Aq, or equivalently 
a. At 4th order in a, there are 2-loop Feynman diagrams that depend logarithmically on 
the ultraviolet cutoff A for the wavenumber of atoms in intermediate states J7|. They give 
an additional momentum-independent contribution to 7^ that is universal. The T-matrix 
element therefore includes the terms 

r 3 = - Bo + 3 ^-^^io 8 ^ + ..., (i2) 

32tt 3 h 6 k K J 

where k is a wavenumber set by the energies of the scattering particles. In the original 
theory, the ultraviolet cutoff is provided by the interaction potential U(r). The ultraviolet 



divergence in (|12|) therefore indicates that there is a momentum-independent contribution 
to 3 — > 3 scattering that depends on the interaction potential. The T-matrix element ([12]) 
can be made independent of the cutoff A if the coefficient Bq depends on A as specified by 
the following renormalization group equation: 

d 3(4tt - 373) m 3 Al 

A dX Bo = ^2^ 7T- (13) 

The parameter Bo could in principle be determined from the numerical solution to the three- 
body problem for the scattering of three atoms interacting pairwise through the potential 
U(r). However, it can equally well be determined from calculations of the ground state 
energy for a large number of atoms interacting through the potential U(r). This is the 
strategy we will follow in this paper. 

Using dimensional analysis, the coefficient of the three-body contact interaction term in 
(§) can be written as 

B = x (144?rc) , (14) 

m 

with c a dimensionless constant. The factor of 1447T has been chosen for later convenience. 
We will refer to c as the three-body contact parameter. 



III. NONUNIVERSAL CORRECTIONS 

In this section, we calculate the leading and next-to-leading nonuniversal corrections to 
the energy density of a homogeneous Bose gas. We also calculate the leading nonuniversal 
corrections to the condensate fraction. 

A. Mean-field corrections to the energy density 

The leading nonuniversal correction to the energy density £ can be obtained by applying 
the mean-field approximation to the Lagrangian £ e g in (H), including the (if)*ifj) 3 term. In 
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the mean-field approximation, ijj is replaced by its ground-state expectation value v, which 
can be taken real- valued. The condensate number density is n = v 2 . The free energy 
density T in the mean- field approximation is simply — £ e fr, with if> replaced by v: 

^nf = ~^v 2 + -A G v A + ^-B v 6 . (15) 
4 3o 

Setting v 2 = n in the energy density £ = T + fin, it becomes 

8x01= -A nl + —B nl + ^(n-no) . (16) 

The chemical potential is obtained from the condition dJ- ' jdv = 0, which guarantees that v 
minimizes the free energy: 

/imf = -xAqUq + Y^B n 2 . (17) 

Substituting the expressions ( |10|) and fllil) for the coefficients into ([TBI) , we have 

2tt/z 2 



fmf = an 

m 



1 + c(167rn a ) + /x (n — n Q ) . (18) 



In the mean-field approximation, the number density n = tp*if) reduces to n = v 2 . Setting 
n = no in (0), the first term gives the familiar universal mean- field energy density: 

2nh 2 2 

Swat = an ■ ( 19 ) 
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The constant c in (|I8D corresponds to a nonuniversal contribution to the constant d under 
the logarithm in the coefficient of na? in (jl]). 

The universal log(na 3 ) term in the coefficient of na 3 in (II) can be obtained through 
complicated two-loop calculations in many-body quantum field theory. However, as pointed 
out by Braaten and Nieto [0] , it follows from the term in fll8|) proportional to c by a simple 
renormalization group argument. If 3 — > 3 scattering amplitudes are to be independent 
of the ultraviolet cutoff A on the wavenumber of atoms in the intermediate states, the 
coefficient c defined by ([H]) must depend logarithmically on the ultraviolet cutoff A in the 
way prescribed by the renormalization group equation (p!3|). The mean-field contribution 
([TED to the energy density depends on the ultraviolet cutoff through the coefficient c(A). 
The total energy density after renormalization must be independent of A. Therefore, there 
must be loop diagrams in many-body perturbation theory that have logarithmic ultraviolet 
divergences that precisely cancel the A-dependence of c(A). The only momentum scale in 
these diagrams is the inverse of the coherence length £ = (167rr7,oa) -1 ^ 2 , which appears in the 
quasiparticle dispersion relation. Thus the logarithmic divergences are of the form log(A£). 
The energy density can be independent of A only if the explicit logarithmic divergences and 
the coefficient c(A) appear in the combination 

c(A) - 47r ~ 3v ^ log(A/Vl67rn a) 
on 

... 4^ — 3^, , . v 47r — 3v^3 , n N / % 

= c (A) log Aa + log(167™ a 3 ) . (20) 

07T D7T 



8 



The three-body contact parameter c need not be defined by a wavenumber cutoff A. For ex- 
ample, one can define a A-independent parameter c by absorbing the logarithm of Aa in (|20"D 
into c(A). Alternatively, one can use dimensional regularization and minimal subtraction to 
define a parameter c(k) that depends on an adjustable renormalization scale k. One could 
even define c by making a specific choice for the universal constant under the logarithm in 
the low-density expansion for the energy per particle in ([J). We will adopt such a definition 
later in this paper. Regardless of the definition, any dependence of a many-body observable 
on c must be accompanied by a logarithm of n a 3 in the combination 

c+ 4 *~ 3 ^ log(167rn a 3 ). (21) 

Replacing c in (^) by this combination and setting n = n, we reproduce the universal 
logarithmic na 3 correction in ([J). 



B. Semiclassical corrections to the energy density 

The next-to-leading order nonuniversal corrections to the energy density can be obtained 
by computing the semiclassical (one-loop) correction to S, including the effects of the effec- 
tive range term and the three-body contact interaction in C e g. The semiclassical correction 
before renormalization is simply the sum of the zero-point energies of the quasiparticles: 

W^w- (22) 

where e(k) is the quasiparticle dispersion relation. This expression is ultraviolet divergent, 
and renormalization is necessary to obtain a finite result. If we insert the Bogoliubov dis- 
persion relation 

h 2 I 

eBog(^) = — kJ k 2 + lQiran , (23) 
2m v 

renormalization can be accomplished by subtracting powers of k from the integrand to make 
it convergent: 
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The three subtractions can be identified with renormalizations of the vacuum energy, the 
chemical potential, and the scattering length, respectively. If we imposed a cutoff A on the 
wavenumber, the three subtractions would be proportional to A 5 , A 3 , and A, respectively. 
The result of evaluating the convergent integral in (^4j) is 

h 2 

£LHY = ^-^(167m a) 5/2 . (25) 

Setting n = n, this reproduces the correction in (P first calculated by Lee, Huang, 

and Yang M . 
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The result (^5|) can also be obtained by using dimensional regularization. The integral 
is generalized from 3 dimensions to 3 — 2e dimensions, evaluated for a value of e for which 
it converges, and then analytically continued to e = 0: 



ft 2 roo i 

£lhy = / dkk 3 - 2e Jk 2 + 167rcm ( 
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(26) 



To be more specific, the integral is separated into two regions: k < k* and k > k*, where k* 
is an arbitrary wavenumber. The contribution from k < k* is evaluated for complex values 
of e satisfying Re e < 2. The contribution from k > k* is evaluated for complex values of e 
satisfying Ree > §. The sum of the two contributions is 



h 2 T(2-e)T f 5 



f "°8A 2r ( -i) 2 ' ('*""■«>*-• < 27 > 

Setting e = 0, we recover (p5|). 

The effective range and three-body contact terms in (^[) change the quasiparticle dis- 
persion relation e(k). To determine e(k), we substitute i[>(r,t) — v + 4>(v,t) into (|3|). The 
quasiparticle Lagrangian Li consists of the terms that are second order in <p. Using the 
mean-field expression (|17|) for the chemical potential and setting v 2 = no, the quasiparticle 
Lagrangian reduces to 
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(167rn a)[l + 3c(167m a 3 )] (Re0) 2 . (28) 
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The quasiparticle dispersion relation e(k) can be determined by Fourier transforming the 
fields Re0 and Im0 and finding the eigenvalues of the resulting 2x2 matrix. The result is 

h 2 



e(k) = —k^J(l + 47m a 2 r s )k 2 + 16vran [l + 3c(167m a 3 )] . (29) 

Setting r s = c = 0, we recover the Bogoliubov dispersion relation in ( 23|) . 

Since the dispersion relation ( |2"9"| ) has the same functional form as the Bogoliubov disper- 
sion relation (^3|) , we can compute the semiclassical energy simply by rescaling the integrand 
and the integration variable k in the dimensionally regularized expression (|26|): 

£ ™ = tJ^^ 1 + ^n a 2 r s y 2 {l67rn Q a[l + 3c(167m a 3 )]} 5/2 . (30) 
Expanding to first order in r s and c, this reduces to 

2nh 2 9 f 32 f— r 16 16 r s l ,„ „ □ sq/ol ,„„ s 

£ sc = an 2 —^16nn a 3 + — c - — - 167m a 3 f' 2 . 31 

m I 157T V L 7T 157T a J J 
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C. Semiclassical corrections to the number density 



The total number density is n = (ifj^ip), while the number density of atoms in the 
condensate is n = \(ip)\ 2 - In the mean-field approximation, the operator ip is replaced by its 
ground-state expectation value v, in which case n = n Q . In the semiclassical approximation, 
we have 



ru 



n + \ 



d 3 k 



Ti z k 2 /(2m) 



+ 



e(k) 



h 2 k 2 /{2 



m) 



(32) 



Inserting the Bogoliubov dispersion relation (|23|) and using dimensional regularization to 
remove the power ultraviolet divergences, this can be written 



^Bos 



n + 
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(33) 



e->0 



Evaluating the integrals and analytically continuing to e = 0, we obtain the classic result of 
Bogoliubov M: 



^0 + 
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24tt 2 



'167m n a) 3/2 



(34) 



In order to compute the complete semiclassical correction to the condensate density, 
including the effects of the effective range and the three-body contact coefficient, we must 
insert the quasiparticle dispersion relation (^) into the integral (^3). Since (|29D has the 
same functional form as the Bogoliubov dispersion relation fl23|) , we can evaluate the integral 
simply by rescaling the integrand and the integration variable k in the expression (j33|): 



n + 



24tt 2 



2(1 + 4:im a 2 r s y 2 - (1 + 47m a 2 r s )~ 1 



(35) 



x 



|l67moa[l + 3c(167moa 3 )] j 
Expanding to first order in r s and c, this reduces to 



3/2 



n ( 1 + — \/lQim a 3 + 
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(36) 



D. Low-density expansion for the condensate fraction 

The semiclassical expression ( j36[) for the number density includes both the leading uni- 
versal corrections and the leading nonuniversal corrections to the mean-field result n = n Q . 
The coefficient c in the (noa 3 ) 3 / 2 correction must be accompanied by a universal logarithmic 
correction in the combination (^). There are universal corrections proportional to (noa 3 ) 
and (noa 3 ) 3 ^ 2 whose coefficients are as yet unknown. Including all corrections through order 
(noct 3 ) 3//2 , the expression for the number density is 
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2tt a 



+ e (167rnoa 



3\3/2 



(37) 



The parameters ce and r s depend on the interaction potential, while d and e are universal 
coefficients. The subscript on ce is a reminder that we will define this parameter by the 
form of the low-density expansion for the energy of the homogeneous gas. 

The low-density expansion for the condensate fraction n^jn can be obtained by inverting 
(|TF|) to obtain no as a function of n. Including all terms through order (na 3 ) 3//2 , the expression 
for the condensate fraction is 



no 
n 



where d! 



1 VlQnria 

3tt 



5 - d'(lQirna 3 



ce + 



67T 



log(167rna ) 



27r a 



+ e' (16vma 3 ) 3 / 2 , 



(38) 



d — „\ and e' 



3tt U ^ 9tt3- 



E. Low-density expansion for the energy density 

The mean-field contribution to the energy density in ( |18D includes the leading nonuni- 
versal correction to the universal mean-field energy density. The semiclassical contribution 
in (^l|) includes the leading universal correction and the next-to-leading nonuniversal correc- 
tions. The coefficient c in the nonuniversal correction must be accompanied by a universal 
logarithmic correction in the combination (|21~D . There is a universal correction proportional 
to n a 3 that was calculated by Braaten and Nieto [I7|. There is also a universal correction 
proportional to (n a 3 ) 3 / 2 that has not been calculated. Including all corrections through 
order (n a 3 ) 3 ^ 2 , the expression for the energy density is 



£ 



-anT, { 1 + 



m 



+fi (n 




167moa 3 + 
4^-3^ 



Ce H ~ log(167rn a ) 



c E + 



log(167rn a 3 ) 



9tt 2 



(167rnoa 3 ) 



16 r s 
157T a 



+ b) (167rn a 3 ) 3/2 



(39) 



The parameters ce and r s depend on the interaction potential, while b is a universal co- 
efficient. The chemical potential u is determined by the condition dT/dv = 0, which is 
equivalent to 



d_ 

dm 



{£ - /in) = 0. 



(40) 



The constant 4/(97r 2 ) in the coefficient of the noa 3 correction in (|39| ) is part of the universal 
coefficient. The rest of that coefficient has been absorbed into the nonuniversal coefficient eg. 
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This represents a specific definition of the three-body contact parameter c whose motivation 
will become evident below in Eq. ([ill). 

The low-density expansion for the energy density can be obtained by using ( |38"D to elim- 
inate n from (|39|) in favor of n. Including all terms through order (na 3 ) 3 / 2 , the expression 
for the energy density is 



1 + 



umf 



32 

lOTT 



+ 



Vl6nna 3 + 



ce H log(lo7rna J 



67T 



(16iina 3 



16 



7T 



4^-3^/3, , „ 
c E H log(167ma ) 

07T 



16 r s 
157T a 



+ 6') (16W) 3 / 2 



(41) 



3tt 97T 3 



where V — b 

We have chosen to define the three-body contact parameter ce so that it includes all the 
universal terms under the logarithm log(167rna 3 ) in the na 3 term in fl41f) . This is a perfectly 
acceptable definition of the three-body contact parameter, because it can be related to any 
other definition by a well-defined calculation. For example, dimensional regularization and 
minimal subtraction can be used to define a three-body coefficient cms( k ) that depends on 

The two-body and three-body coupling constants 
and Bq up to factors of h. The energy density 



an arbitrary renormalization scale k |L7 
g and g% (k) in Ref. 



17H are equal to A 

was calculated to order na 3 using dimensional regularization and minimal subtraction in 
Ref. ||17| . By equating the coefficient of the (167ma 3 ) term in the expression (fy]) with the 
corresponding coefficient in Ref. |T^], we find 



c E + 



4tt-3\/3 

67T 



log(167rna 3 



cms («) + 



4tt-3\/3 

67T 



16rma 
log — + 0.80 



(42) 



Thus the parameter ce defined by the expression (|4l]) is equal to the running parameter 
cms( k ) at a scale of order 1/a: 



ce = c ms {k = 1.5/a) 



(43) 



IV. DIFFUSION MONTE CARLO RESULTS 

In this section, we present the diffusion Monte Carlo results of Giorgini, Boronat, and 
Casulleras for the energy and the condensate fraction of the homogeneous Bose gas. We 
also compute the effective ranges for each of the model potentials they considered. 



A. Model potentials 

In the low-density limit, a homogeneous Bose gas exhibits universal behavior that is 
determined completely by the S-wave scattering length a. However, at any nonzero density, 
its properties will also depend at some level on aspects of the interactions between atoms 
other than the scattering length. These nonuniversal effects scale as higher powers of n and 
become negligible in the low-density limit. 
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If Ref. U, the authors considered 4 different two-body potentials as models for the 
interactions between the bosons. These potentials all have the same scattering length a. 
The potentials are 

• the hard-sphere potential (HS): 

HS )(r) = oo, r<a, 

= 0, r > a. (44) 

• two soft-sphere potentials (SS-5 and SS-10): 

y (ss) (r) = V , r<R, 

= 0, r>R. (45) 

The radii are R = 5 a for SS-5 and R = 10 a for SS-10. The heights of the two potentials 
are V = 0.031543 h 2 /(ma 2 ) for SS-5 and V = 0.0034084 h 2 /{ma 2 ) for SS-10. 

• a hard-core square- well potential (HCSW): 

= -V , R c <r<R, 

= 0, r>i2. (46) 

The outer radii of the hard core and the square well are R c = a/50 and R = a/10, 
respectively. The depth of the square well is V = 412.815 h 2 /(ma 2 ). The potential 
supports a single two-body bound state with energy -1.13249 h 2 /(ma 2 ). 

Nonuniversal effects in the Bose gas depend on aspects of the interactions between the 
bosons other than the S-wave scattering length. One of the most important is the effective 
range for S-wave scattering r s , which is defined by the low-momentum expansion (§) of the 
phase shift for atom-atom scattering. The effective range for the hard-sphere potential is 
r s = 2a/3. For the soft-sphere potentials, the S-wave scattering length is 

tanh(fc i?) 

a = R , (47) 

fc 



where A; = ymVo/h 2 . The effective range is 

R 3 1 

Its numerical value is r s = —4.96372 a for SS-5 and r s = —29.936 a for SS-10. For the HCSW 
potential, the S-wave scattering length is 

a = R _ ^MR-R c )} _ (49) 

Kq 



The effective range is 



„ R 3 R c (R-a) 2 R c -a 
3a 2 a 2 kia 2 



Its numerical value is r, = 0.113493 a. 
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SS-10 


HCSW 


i x icr 6 


1 00427(80) 


1 00427(80) 

i . U Ut:*£j t \ UU 1 


1 00427(80) 


1 00427(80) 


5 x icr 6 


1 00952(16) 








8 x icr 6 


1 01262(99) 




1 01163(99) 




1 x 10~ 5 


1 01382(80) 


1 01382(80) 


1 01302(80) 


1 0162(16) 


1.2 x icr 5 


1 01528(66) 




1 01395(66) 




5 x 10~ 5 


1.02957(48) 








i x icr 4 


1 04326(80) 


1 04167(80) 


1 03689(80) 


1 04565(80) 

XaU^UUUlUU j 


r x in -4 


1 09499(64) 








1 x 10~ 3 


1 1332(16) 


1 11011(80) 


1 07907(80) 


1 1380(40) 


5 x ID -3 

u A J.U 


1 29775(95) 








i x i(r 2 


1.42921(80) 








5 x 10~ 2 


2.1295(16) 








1 x KT 1 


2.8855(56) 








0.166 


3.9597(96) 








0.244 


5.446(33) 









TABLE I. Energy normalized to the universal mean-field energy for the HS, SS-5, SS-10, and 
HCSW potentials from the diffusion Monte Carlo calculations of Ref. ||. The statistical errors are 
given in parentheses. 

B. Energy and Condensate Fraction 

In Ref. , the authors used the diffusion Monte Carlo method to calculate the properties 
of a homogeneous gas of bosons interacting through each of the four model potentials. In 
the case of the potentials HS, SS-5, and SS-10, the homogeneous gas is the ground state of 
the system. In the case of the potential HCSW, the homogeneous gas is a metastable state 
in which there are no two-body bound states. 

They calculated the energy per particle E/N and the condensate fraction N /N as 
functions of na 3 , with na 3 ranging over several orders of mag nitude from 1(T 6 to 0.244. 
Their results are reproduced in Tables | and |J. The condensate fraction was not calculated 
for the HCSW potential. Most of the values in Table | were given in a table in Ref. ||. 
The energy values for na 3 = 8 x 10~ 6 and na 3 = 1.2 x 10 -5 and all the condensate values 
in Table [TI] were provided to us by the authors. For convenience, we have normalized E/N 
to the universal mean-field result E um f/N = (2-kTi 2 /m)an. The errors in parentheses are the 
statistical errors from the diffusion Monte Carlo. The number of particles in the simulation 
box was N = 500, which the authors claimed is large enough that finite-size effects are well 
below the statistical errors. 

V. ANALYSIS 

In this section, we obtain crude determinations of the three-body contact parameter ce 
defined by fl4"T| ) for each of the interaction potentials by analyzing the diffusion Monte Carlo 
data for the energy. We also use the diffusion Monte Carlo data for the HS potential to 
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0.998(1) 


0.997(1) 


0.997(2) 


5 x 1(T 6 


0.997(1) 






1 x 1(T 5 


0.996(2) 


0.997(2) 


0.996(1) 


5 x 1(T 5 


0.991(2) 






1 x 1(T 4 


0.987(3) 


0.989(2) 


0.992(1) 


5 x 10~ 4 


0.967(3) 






1 x 1(T 3 


0.948(5) 


0.971(3) 


0.986(1) 


5 x 10~ 3 


0.876(6) 






1 x 1(T 2 


0.803(3) 






5 x 1(T 2 


0.501(5) 






1 x 


0.277(3) 






0.166 


0.109(4) 






0.244 


0.031(4) 







TABLE II. Condensate fractions for the HS, SS-5, and SS-10 potentials from the diffusion 
Monte Carlo calculations of Ref. ||. The statistical errors are given in parentheses. 



obtain a crude determination of the universal coefficient d' in the low-density expansion ( |38| ) 
for the condensate fraction. 



A. Three-body contact parameters 



In Section |V|, we have calculated the effective range r^ pot ) for each of the four model 
potentials used in Ref. 0. As a consequence, the only unknowns in the expression (|4i~D for 
the energy density are the three-body coefficient c^ ot \ which depends on the interaction 
potential, and the universal coefficient V . At small densities, the contribution of V is sup- 
pressed by a power of \/na 3 compared to the contribution of Cg 0t \ We define a function 



C(p°t)( 



na" 



that in the limit na 3 



approaches c^ >ot ' ) via 



C ipot \na 3 ) 



1 



16nna' 3 



£(p°t) 



umf 



32 

loTT 



Vl6 



una — 



4tt-3\/3 

67T 



(167ma 3 ) log(167ma 3 



• (51) 



The approach to the limiting value c (pot) as 



(J (pot) 



„(pot) 



+ 



16 



7T 



(pot) . 47T - 3a/3 3 

Cr H log(167rna J ) 

67T 



na 3 —> is given by 

Ig r (pot) 



157r a 



+ b'\ (167rna 3 ) 1/2 . (52) 



For sufficiently small values of na 3 , the approach of the function C^ pot ^ to the constant value 
C£ >ot ' ) should be dominated by the logarithmic term in (|52|). 

We first consider the HS potential. In Fig. [1], we show the values of C^ 3 ' calculated 
from the diffusion Monte Carlo data in Table |I[ The curves give some idea of the expected 
size of the corrections from next order in Vna 3 . The dashed and dot-dashed curves show 
the corrections from the logarithmic term and the effective range term, respectively. The 
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FIG. 1. The value of C^ HS ^ for the data points from Table |[ The dot-dashed curve shows the 
corrections from the effective range term. The dashed curve in shows the corrections from the 
logarithmic term at the next order. 



correction from the b' term in (52) has the same density dependence as the effective range 
term, but its normalization in unknown. The logarithmic correction will be the largest in 
the region na 3 < 1.2 x 10~ 5 unless \b'\ > 14.8 or 



,(HS), 



> 2.9. For na 3 < 1.2 • 10" 



the 



statistical errors in the data are larger than the expected size of the corrections from higher 
orders. We can therefore obtain a rough determination of c^ 1 ^ by averaging C^ 1 ^ at the 5 
density values na 3 < 1.2 x 10~ 5 weighted by the inverse-squares of the statistical errors. The 
lowest data point could equally well be omitted because of its large error bar. The resulting 
value is 



„(HS) 



-0.9 ±0.5. 



(53) 



The error was estimated by adding the inverses of the statistical errors in quadrature and 
then taking the reciprocal. The value fl53p is dominated by the data point at na 3 = 5 x 10~ 6 , 



because it has such a small statistical error. The systematic error from the \fncfi corrections 



in (52) is expected to be smaller than the magnitude of the logarithmic correction in ([521) , 
which is 0.26 at na 3 = 5 x 10 -6 . 

We also performed a global x 2_ fit of the data using ([HP through order na 3 . We assigned 
a theoretical error ±(167ma 3 ) 3 / 2 to quantify the error from leaving out the unknown higher 
order terms in the expansion. We computed the x 2 f° r the difference between the data 
in Table | and the theoretical expression ( f41~D through order na 3 using the error obtained 
by adding the theoretical error and the statistical error from Table [I] in quadrature. The 

/TTO\ 

preferred value is c E = 0.8. In this case, the data points at higher density favor a positive 
value for c^ S ^ in spite of their larger theoretical error. 
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We next consider the SS-10 potential, for which there are six data points in Table j} In 
Fig. |2|, we show the values of C^ ss ~ 10 ^ calculated from the diffusion Monte Carlo data for the 
SS-10 potential in Table [I]. The curves are the same as in Fig. [I], except that the effective 
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FIG. 2. The values of C^ ss 10 ) for the data points from Table |[ The curves are as in Fig. H. 

range correction is much larger, because 7^ ss ~ 10 ) is larger than rf 1 ^ by about a factor of 45. 
For na 3 < 1.2 • 10~ 5 , the statistical errors in the data are larger than the expected size of the 
corrections from higher orders. We can therefore obtain a rough determination of c( ss-1 °) by 
averaging C'( ss_10 ) at the 4 density values na 3 < 1.2 x 1CF 5 weighted by the inverse-squares 
of the statistical errors. The resulting value is 



JSS-10) 



■1.6 ±0.8. (54) 



For both the SS-5 and HCSW potential, there are only four data points in Table |. In 
Fig. ^ we show the corresponding values of (7^ ss_5 ^ and (7( HCSW ). The expected error from 
higher orders is smaller than the statistical error for the 2 data points with na 3 < 10 -5 . If 
we use the value of C^ pot ^ at na 3 = 10~ 5 as an estimate for Cg 0t \ we obtain c^ S 5 ^ = 0.2± 1.6 
and 4 HCSW) = 4.9 ± 3.2. 

We have given a rough determinations of the three-body contact parameter c^ ot ^ for 
each of the potentials HS, SS-10, SS-5, and HCSW. The error bars are large enough that we 
cannot identify a deviation from universality in the three-body coefficients. More accurate 
determinations would require higher statistics in the diffusion Monte Carlo calculations. 
It might also require increasing the number of particles in the simulation box in order to 
keep finite-size effects smaller than the statistical errors. Ideally, one would like to have 
high enough statistics at several values of na 3 between 10 -4 and 10~ 6 so that the rate of 
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FIG. 3. The values of C( ss ~ 5 ) (left plot) and C( HCSW ) (right plot) for the data points from 
Table |. The curves are as in Fig. [I]. 



approach of C^ pot \na 3 ) to the constant value c^ >ot ' ) could also be used to determine the 
universal coefficient b' in (1521). 



B. Universal 2nd-order correction to the condensate fraction 



In the expression (j38j) for the condensate fraction, the only unknowns are the three-body 



coefficient c^ ^ and the universal coefficients d! and e'. We define a function D^ pot h 
that approaches the universal coefficient d' as the density approaches zero by 



£)(P ot )| 



na 



n J 



(pot) 



1 H v IQuna 3 

3tt 



167rna 3 

The approach to the limiting value d! as na 3 — > is given by 

]_ r (p°t) 



(55) 



£)( P Ot) 




c (pot ) + 4 7 r-3v / 3 log(16W) 

67T 



27T a 



+ e' (16vma 3 ) 1/2 . (56) 



In Fig. ^, we show the values of _D( HS ) calculated from the diffusion Monte Carlo data for 
the HS potential in Table [T|. The curves give some idea of the expected size of the corrections 
from next order in Vna 3 . The dashed and dot-dashed curves show the corrections from the 
logarithmic term and the effective range term, respectively. The correction from the e! term 
in (|56|) has the same density dependence as the effective range term, but its normalization in 
unknown. The logarithmic correction should be the largest in the region na 3 < 10~ 4 , unless 
\e'\ > 2.0 or \cf S) \ > 2.1. For na 3 < 10 4 , the statistical errors in the data are larger than 
the expected size of the corrections from higher orders. We can therefore obtain a rough 
determination of d' by averaging _D( HS ) at the 5 density values na 3 < 10~ 4 weighted by the 
inverse-squares of the statistical errors. The resulting value is 
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FIG. 4. The values of Z)( HS ) calculated from the diffusion Monte Carlo data in Table |Il|. The 
dashed line shows the corrections from the logarithmic term while the dash-dotted line shows the 
corrections from the effective range. 



d' = -0.5±0.5. (57) 

The systematic error from the yna 3 corrections in (j56|) is expected to be smaller than the 
size of the logarithmic correction in (^), which is 0.14. We have also determined d' by 
using a global x 2_n ^ to the data of Eq. fl3~7|) through order na 3 , using a theoretical error 
of ±(16nna 3 ) 3 / 2 that is added in quadrature to the statistical error. The preferred value is 
d! — 0, which agrees with Eq. (|57D to within the error bars. 

For both the SS-5 and SS-10 potentials, there are only 4 data points for the condensate 
fraction in Table 0. The condensate fraction was not calculated in Ref. for the HCSW 
potential. In Fig. we show the values of _D( SS ~ 5 ) and _D( SS-1 °) calculated from the diffusion 
Monte Carlo data in Table O. The curves are the same as in Fig. fj, except for the different 
value of the effective range. The expected error from higher orders is smaller than the 
statistical error for the 2 data points with na 3 < 10 -5 . If we use the value of £> (pot) at 
na 3 = 10~ 5 as an estimate for d f , we obtain d' = —3.5 ± 4.0 from the SS-5 potential and 
d! = — 1.5 ± 2.0 from the SS-10 potential. These values are consistent within errors with the 
value ( |5TD obtained from the condensate fraction for the hard-shere potential. 



VI. SUMMARY 



Effective field theory is an ideal tool for the theoretical study of nonuniversal effects in the 
homogeneous Bose gas. The central idea of effective field theory is to exploit the separation 
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FIG. 5. The values of Z?( ss_5 ) (left plot) and Z?( ss ~ 10 ) (right plot) calculated from the diffusion 
Monte Carlo data in Table |IJ The curves are as in Fig. The first data point at na 
29.8 ± 19.9 for the SS-5 potential and 29.8 ± 39.8 for the SS-10 potential. 



3 - 1(T 6 is 



of scales in a physical system. In a low-density Bose gas, the important length scales include 
the interparticle spacing, which is proportional to n" 1 ^ 3 , and the coherence length, which is 
proportional to (na)" 1 / 2 . These are both much larger than the physical size of the atoms, and 
in typical experiments they are also much larger than the scattering length a. The number 
of atoms per cubic scattering length na 3 is therefore a small parameter, and it allows for 
a controlled, perturbative expansion of low-density observables in yna 3 . This systematic 
expansion is conveniently implemented using effective field theory. In low-energy observables, 
all the sensitivity to the details of the interactions between atoms at short distances is 
captured in low-energy constants, such as the scattering length, the effective range, and the 
effective three-body contact interaction. Universal results that depend only on the scattering 
length can also be obtained using more traditional many-body methods, but the calculations 
can be simplified using effective field theory. However, the real power of effective field theory 
becomes apparent in the systematic calculation of nonuniversal corrections. 

Effective field theory can be used to easily determine at what order in the low-density 
expansion various nonuniversal effects will enter. The scattering length is associated with 
the (ip*ip) 2 term in the local effective Lagrangian (0). The low-energy parameters associated 
with nonuniversal effects can be identified with terms in the effective Lagrangian that are 
higher order in ip*ifj or in V. In general, there are infinitely many such terms, but only a 
finite number contribute to any observable at a given order in the low-density expansion. 
The order in y/n at which a given term contributes is determined by the powers of ip*ip 
and V and by the order in the loop expansion at which it first enters. For example, if we 
consider the energy density, every power of ip*ij) contributes a factor of n. Every power 
of V contributes a factor of \fna. Finally there is a quantum suppression factor of Vna 3 
associated with each order in the loop expansion. For example, the three-body contact 
term (ip*ip) 3 gives a mean-field contribution proportional to n 3 . The effective range term 
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\S7(ip*ip)] 2 does not contribute in the mean-field approximation, because the gradients vanish 
for a homogeneous gas. Its first contribution is at the semiclassical (one-loop) order and is 



proportional to ^/na^{^/na) 2 n 2 ~ n 7 / 2 . Thus the leading nonuniversal contribution comes 
from the three-body contact interaction. The effective range enters at one order higher in 
y/n. At order n A , there is also a nonuniversal contribution to the energy density associated 
with the 4 term in the effective Lagrangian. 

We have calculated the leading and next-to-leading nonuniversal corrections to the energy 
density and the leading nonuniversal corrections to the condensate fraction for a homoge- 
neous Bose gas. These nonuniversal corrections are determined by the effective range r s and 
the three-body contact parameter c. These parameters are in general different for different 
potentials, even if they have the same scattering length a. For a given potential, the effective 
range can be extracted from the amplitude for 2-^2 scattering in the standard way. The 
parameter c could in principle be determined from a three-body observable, such as the 
amplitude for 3^3 scattering in the vacuum. However, it can equally well be obtained 
from a low-density observable. In fact, we chose to define the parameter ce using the form 
of the low-density expansion for the energy density in PT|). 

We analyzed the nonuniversal effects in the diffusion Monte Carlo results of Giorgini, 
Boronat, and Casulleras || for a homogeneous Bose gas with four different interaction po- 
tentials. We attempted to extract the unknown three-body coefficient eg from the diffusion 
Monte Carlo data for the energy density. We found that the statistical errors at small values 
of na 3 allowed only a crude determination of this parameter. Our values for the hard-sphere 
and SS-10 potentials are given in (|53| ) and fl54|), respectively. Within the error bars, we were 
unable to see a deviation from universality in the three-body coefficients. To determine 
ce precisely enough to see the nonuniversality, one would need higher statistics at several 
densities in the range 10~ 4 < na 3 < 10~ 6 . 

We also attempted to extract from the diffusion Monte Carlo data the universal coefficient 
d' in the low-density expansion ( |3"8| ) for the condensate fraction. The value obtained from 
the data for the hard-shere potential is given in (f)7]). Again we found that the statistical 
errors at small values of na 3 allowed only a crude determination of this coefficient. In order 
to get a useful determination of d', one would need higher statistics at several densities in the 
range < na 3 < 10~ 5 . The coefficient d' can also be calculated from Feynman diagrams 
using the methods used in Ref. (I7j to compute the na 3 corrections to the energy density. 
Such a calculation would provide a useful check on the accuracy of diffusion Monte Carlo 
results. 

The achievement of Bose-Einstein condensates in atomic vapors has opened up the pos- 
sibility for experimental study of Bose gases. As the experiments become more and more 
precise, it will be important to understand nonuniversal effects. Effective field theory pro- 
vides a perturbative framework within which nonuniversal effects can be systematically 
calculated in the low density limit in terms of a few low-energy parameters. These effects 
can also be calculated nonperturbatively using diffusion Monte Carlo simulations. More 
extensive diffusion Monte Carlo calculations would be extremely valuable for testing the 
understanding of nonuniversal effects provided by effective field theory. 
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